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ABSTRACT 

For the first time, we investigate the X-ray/infrared (IR) correlation for star-forming 
galaxies at z ~ 1, using SPIRE submm data from the recently- launched Herschel Space 
Observatory and deep X-ray data from the 2 Ms Chandra deep field north (CDFN) sur- 
vey. We examine the X-ray/IR correlation in the soft X-ray (SX, 0.5-2 keV) and hard 
X-ray (HX, 2-10 keV) bands by comparing our z ~ 1 SPIRE-detected star- forming 
galaxies (SFGs) to equivalently IR- luminous (Lir > 10 10 L ) samples in the local/low 
redshift Universe. Our results suggest that the X-ray/IR properties of the SPIRE 
SFGs are on average similar to those of their local counterparts, as we find no evi- 
dence for evolution in the isx/^m and Ljjx/iiR ratios with redshift. We note how- 
ever, that at all redshifts, both £gx/£m and i'Hx/i'iR are strongly dependent on IR 
luminosity, with luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs, 
Ltr > 10 11 Lq) having up to an order of magnitude lower values than normal infrared 
galaxies (Lir < 10 11 L Q ). We derive a Lsx-^ir relation and confirm the applicability 
of an existing Lhx-^ir relation for both local and distant LIRGs and ULIRGs, con- 
sistent with a scenario where X-ray luminosity is correlated with the star-formation 

rate fST?m 
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1 INTRODUCTION 

For at least two decades, the infrared/X-ray connection has 
been an important point of focus in extragalactic surveys. 
The various X-ray missions (Einstein, ROSAT, ASCA, Bep- 
poSAX, Chandra and XMM-Newton) together with the in- 
frared observatories (IRAS, ISO and Spitzer) have estab- 
lished the X-ray/infrared synergy not only as a tool for 
studying the active galactic nucleus/host galaxy interplay, 
but also star-formation and stellar evolution (a few examples 
of the multitude of papers on this topic: Eales & Arnaud 
1988; Fabbiano 1988; Fabbiano, Gioia & Trinchieri 1988; 
Rieke 1988; Green et al. 1989; Rephaeli et al. 1991; Green 
et al. 1992; Boiler et al. 1992; Barcons et al. 1995; Rephaeli, 
Gruber & Persic 1995; Imanishi & Ueno 1999; Severgnini et 
al. 2000; Risaliti et al. 2000; Cavaliere, Giacconi & Menci 
2000; Fadda et al. 2002; Gandhi & Fabian 2003; Manners et 
al. 2004; Panessa et al. 2005; Netzer et al. 2005; Franceschini 
et al. 2005; Hickox et al. 2007; Treister et al. 2009; Trichas 
et al. 2009; Lutz et al. 2010; Park et al. 2010; Goulding et 
al. 2011). 

Emission in the infrared (IR) to millimetre (~ 8- 
1000 /an) part of the spectrum is attributed to galactic dust 
being heated to a range of temperatures by incident UV and 
optical photons. In cases where an active galactic nucleus 
(AGN) is present, emission from the torus (or other axisym- 
metric dust/gas configurations, e.g. Efstathiou & Rowan- 
Robinson 1995) can also play a prominent role at short in- 
frared wavelengths (<20/xm, rest-frame), as torus dust is 
heated to near-sublimation temperatures by the AGN radi- 
ation field. On the other hand, flux longwards of ~20/im 
is strongly correlated with the star formation rate (SFR), 
although in the case of quiescently star-forming galaxies, 
emission from dust heated by older stellar populations might 
also contribute significantly to the far-IR (e.g. Bendo et al. 
2010, Calzetti et al. 2010). Nevertheless, in galaxies with 
high total infrared luminosities (Lir >10 10 L0) where large 
dust masses are implied, (i) the source of heating is intense 
and (ii) the dust replenishment is continuous, making IR 
emission an excellent tracer of young star-forming/starburst 
regions. Moreover, contamination from an AGN (if present) 
is minimal in the far-IR/submm part of the spectrum. Re- 
cently, Hatziminaoglou et al. (2010) showed that the submm 
colours of AGN hosts are indistinguishable from those of the 
general submm-detected population. 

Emission in the high energy part of the spectrum is a 
combination of both AGN (if present) and star-formation 
related processes. AGN are the most powerful X-ray emit- 
ters, with X-rays originating in the hot corona surrounding 
the accretion disk through inverse Compton scattering of 
softer energy photons (Sunyaev & Titarchuk 1980), gener- 
ally amounting to a few per cent of the AGN energy budget 
(Elvis 1994). On the other hand, the applicability of X-rays 
as a star formation tracer is linked to a galaxy's star forma- 
tion history and interstellar medium (ISM) conditions (e.g. 
Ghosh & White 2001; Fabbiano et al. 2004). Star- forming 
galaxies have complex 0.5-10 keV spectra, with contribu- 
tions from a low temperature soft component (kT <1 keV) 
and a harder component which is either predominantly 
thermal (kT ~5-10 keV) or non-thermal (power-law index 
r~1.5-2) or a combination of both (Persic & Rephaeli 
2002). The soft component is primarily thermal emission 



from gas in the ISM heated to X-ray temperatures by stel- 
lar winds and supernovae, with a contribution from point 
sources of the order of 30 per cent, according to Chandra ob- 
servations of local galaxies (e.g. Fabbiano, Zezas & Murray 
2001; Zezas, Ward & Murray 2003). On the other hand, more 
than 60 per cent of the hard component is resolved into point 
sources and associated with X-ray binaries (e.g. Griffiths et 
al. 2000). High mass X-ray binary (HMXB) spectra can be 
represented by a power law up to about 20keV, whereas 
the spectra of low mass X-ray binaries (LMXB) are a com- 
bination of thermal bremmstrahlung and black-body emis- 
sion in the <20keV region (e.g. Christian & Swank 1997, 
see also Persic & Rephaeli 2002). HMXBs, where the com- 
panion star is massive and short-lived, are direct tracers of 
recent star-formation, whereas LMXBs have lifetimes of the 
order of 10 Gyr and are hence more appropriate tracers of 
stellar mass (e.g. Ptak et al. 2001; Grimm, Gilfanov & Sun- 
yaev 2002, 2003). Although the relative timescales of X-ray 
emitting processes play an important role in how one inter- 
prets the star-formation history of quiescently star-forming 
galaxies or starburst (SB) systems, X-rays are a good over- 
all tracer of the first ~ 30 Myr of star-forming activity (e.g. 
Mass- Hesse, Oti-Floranes & Cervino 2008), with the rela- 
tion between star formation rate (SFR) and full band X-ray 
luminosity found to be linear in the >1 Mq yr _1 SFR regime 
(Gilfanov, Grimm & Sunyaev 2004). 

As both X-ray and IR emission trace star formation, 
one would expect these quantities to be correlated in the 
absence of an AGN. The X-ray/IR correlation, described 
as log (Lx)=A + B log (Lir) , has been studied extensively in 
both the soft (0.5-2 keV) and hard (2-10 keV) energy bands 
with samples of local/low-redshift galaxies (eg. Griffiths & 
Padovani 1990; David, Jones & Forman 1992; Ranalli, Co- 
mastri & Setti 2003; Franceschini et al. 2003; Rosa-Gonzalez 
et al. 2007). These studies have calculated different values 
for the slope of the soft X-ray/IR relation with B roughly 
ranging from 0.7 to 1. A linear relation (B=l) is a strong 
indication that X-ray and IR emission directly trace the star 
formation rate, however deviations from linearity do not nec- 
essarily suggest the contrary. The large range in computed 
slopes could be a consequence of the small number statistics 
in the samples studied, in combination with the large scatter 
in the galaxies' physical properties, such as the efficiency of 
energy coupling to the ISM which heats the X-ray emitting 
gas, varying contributions from point sources, as well as the 
amount of absorption undergone by the soft X-rays due to 
gas/dust columns in the line of sight. In that respect, hard 
X-rays are a more reliable tracer of star- formation, as the 
ISM is mostly transparent at energies above ~ 2keV, except 
for dense, massive molecular cloud environments (Grimm, 
Gilfanov & Sunyaev 2002; Ranalli, Comastri & Setti 2003; 
Persic et al. 2004). Indeed, hard X-ray emission linearly 
traces the star formation rate (e.g. Grimm, Gilfanov & Sun- 
yaev 2003; Persic et al. 2004), although because of the lower 
sensitivity of detectors in the hard X-rays, this band is still 
not as extensively used. In addition, whereas soft X-rays 
suffer little contamination from X-ray binary emission, to 
use hard X-rays successfully, the relative contribution of 
HMXBs and LMXBs needs to be quantified. This is easier 
for SB environments where HMXBs dominate but for more 
quiescently star-forming systems, LMXBs, which trace the 
stellar mass, can contribute substantially, introducing non- 
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linearities (e.g. Grimm, Gilfanov & Sunyaev 2002; Persic et 
al 2004). 

So, what is the role of the X-ray/IR correlation in ex- 
tragalactic surveys? Apart from the obvious motivation, i.e. 
examining the physical properties of different galaxy popula- 
tions, one of its key applications could be to discriminate be- 
tween AGN-dominated and star-formation dominated sys- 
tems. The ability to discriminate between the two is par- 
ticularly significant in X-ray surveys, where AGN dominate 
the bright end of the number counts, but the star-forming 
galaxy population is expected to outnumber AGN at faint 
fluxes (fo.5-2kov ~ 10 -17 ergs -1 cm -2 ; e.g. Hornschemeier et 
al. 2002, 2003; Bauer et al. 2004; Brandt & Hasinger 2005; 
see also Georgakakis et al. 2003; 2004; 2006ab; 2007; 2008). 
As previously shown (e.g. Georgakakis et al. 2007), the X- 
ray/IR correlation is a useful way of separating starbursts 
from AGN. However, up to now, such studies have used rela- 
tions applicable to normal IR galaxies in the local Universe. 
What remains to be settled is whether there is a change 
in the characteristics of the X-ray/IR relation with redshift 
and/or SFR. As noted earlier, a linear relationship for star- 
forming galaxies is a strong indication that X-ray and IR 
emission directly trace the star formation rate. As a result, 
determining the X-ray/IR relation would play an important 
role in calibrating X-ray emission as a star formation in- 
dicator, providing an alternative census of the global SFR 
density. 

Although various extragalactic X-ray studies have tar- 
getted cosmologically significant redshifts (e.g. Norman et 
al. 2004; Ptak et al. 2007; Georgakakis et al 2007; Lehmer 
et al. 2008), the X-ray/IR correlation has been examined 
only for local galaxies (e.g. Ranalli, Comastri & Setti et 
al. 2003) and X-ray selected samples at moderate redshifts 
(Grimm, Gilvanov & Sunyaev 2002; Rosa-Gonzalez et al. 
2007), mainly due to sensitivity limitations in X-ray sur- 
veys. For the first time, we are able to investigate the X- 
ray/IR correlation at z~l, for submm-bright sources with 
high star-formation rates, using deep X-ray data from the 
Chandra 2 Ms CDFN survey and sub-mm data from SPIRE 
(Griffin et al. 2010), the sub-mm bolometer array on the 
recently-launched Herschel Space Observatory^] (Pilbratt et 
al. 2010). Our aim is to compare the X-ray properties of 
z ~1 IR-luminous galaxies to their counterparts in the local 
(z < 0.1), moderate (z ~ 0.6) and high (z ~ 2) redshift 
Universe. The paper is laid out as follows: in Section [2] we 
present the sample and in Sections [3] and [4] we examine the 
soft and hard X-ray/IR relations. Our summary and conclu- 
sions are presented in Section [5] Throughout this paper, we 
adopt a concordance consmology of Ho=70kms _1 Mpc _1 , 

n M =i-n A =o.3. 



2 THE DATA AND SAMPLE SELECTION 

2.1 Optical and Infrared data 

The infrared data come from the Herschel /SPIRE Science 
Demonstration Phase (SDP) observations of the Great Ob- 



1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 



servatories Origins Deep Survey (GOODS)-North, taken as 
part of the Herschel Multi-tiered Extragalactic Survey (Her- 
MES) guaranteed-time key programmqj (Oliver et al. in 
prep.). We focus on the region of overlap (0.08 deg 2 ) be- 
tween the HerMES GOODS-N coverage and the deep 24 fj,m 
Spitzer /MIPS (Rieke et al. 2004) observations (Magnelli et 
al. 2009). Source extraction on Herschel /SPIKE 250, 350 
and 500 /im image^j is performed on the IRAC-3.6/im po- 
sitions of the 5 a 24 fj,m GOODS-N sources as described in 
Roseboom et al. (2010). We find 84 sources with /250yjm 
above the 3<r confusion of 17.4 mjy (see Nguyen et al. 
2010 for the SPIRE confusion limits), hereafter referred to 
as the full SPIRE sample. This method of source extrac- 
tion on already known positions is widely used and enables 
identifications of secure counterparts over the whole spec- 
tral energy distribution (SED). In this case however, its 
significant advantage, lies in its ability to effectively deal 
with source blending in the SPIRE bands. Due to the large 
SPIRE beams of 18.1, 24.9 and 36.6 arcsec (FWHM) at 250, 
350 and 500 ^m respectively (Nguyen et al. 2010), many 
SPIRE detections are composed of confused sources. By us- 
ing IRAC-3.6/^m positions of 24/^m sources, we are able to 
assign SPIRE fluxes to each 24/im source in the blended 
composit and hence obtain the aforementioned 'clean' full 
SPIRE sample of 84 sources. A detailed report on the re- 
liability and completeness of this method is presented in 
Roseboom et al. (2010), where it is judged in the context 
of results from the sub-mm number counts and cosmic in- 
frared background. Roseboom et al. demonstrate that for 
deep 24/im surveys and bright SPIRE flux densities, the 
250 /mi-detected, 24 /im sources are essentially the objects 
that make up the 250 fim population. 

Following the same method of source extraction as 
for the SPIRE bands we also obtain photometry at 
Spitzer /MIPS 70 and 160 /im in order to better sample the 
infrared SED. We also retrieve R-band magnitudes, by cross- 
matching the IRAC-3.6/im positions of the full SPIRE sam- 
ple with the Capak et al. (2004) GOODS-N optical cata- 
logue, within a radius of 1.5 arcsec. The association rate is 
90 per cent. 



2.2 X-ray data 

The X-ray observations are from the 2 Ms Chandra Deep 
Field North (CDFN) survey. The methodology for data re- 
duction, source detection and photometry estimates is de- 
scribed in Laird et al. (2009) and the final data products 
are made available by the Imperial College (IC) teanjj. The 
IC catalogue consists of X-ray sources with a Poisson signif- 
icance of < 4 x 10~ 6 (equivalent to > 4.5a in the case of a 
normal distribution) detected independently in four energy 
bands, full (0.5-7keV), soft (0.5-2 keV), hard (2-7keV) and 
ultra-hard (5-7 keV). Soft and hard X-ray fluxes are esti- 
mated in the 0.5-2 and 2-10 keV energy intervals, adopting 



2 hermes.sussex.ac.uk 

3 The data presented in this paper will be released 
through the Herschel Database in Marseille HeDaM 
(hedam.oamp.fr/HerMES) 

4 http://astro.ic.ac.uk/research/data-products-chandra-surveys 
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Table 1. Table of stacking results for the X-ray non-detections in the redshift sub-sample of 67 sources, grouped into LIRGs and ULIRGs. 
Columns 2 and 3 show the mean total infrared luminosity (with 1 a uncertainties) and mean redshift of the group. Columns 4-7 display 
information for the soft and hard band respectively, /gx and /hx (ergs -1 cm -2 ) are the soft and hard stacked fluxes of each group. 
The significance of each stacked signal is shown in brackets, where a refers to the background level. Corresponding luminosities Lsx an< l 
Lux (ergs - 1 ) are estimated using the mean redshift of the group. Note that for the ULIRGs, there is no stacked signal and therefore 
we quote the 3ct upper limit. 



Group 



Soft band: 0.5-2 keV 
log(L m ) (2) / sx logLsx 



fa 



Hard band: 2-10 keV 

log L H x 



LIRGs (18 sources) 11.55 



+0.23 
0.51 



ULIRGs (5 sources) 12.341°-^ 



0.62 


2.37(±0.45)xl0~ 17 

(9.3 cr) 


40.56lH 8 9 


1.34 


2.43(±0.74) xl0~ 17 
(6 a) 


4i.38±g:iS 



7.08(±2.5)xl0 -17 41.04 

(4.3 a) 

<8.03 xlO -17 <41.9 

(<3<r) 



+0.13 
0.2 



a photon index of F — 1.4 and using the Bayesian methodol- 
ogy described in Laird et al. (2009). We cross- match the IC 
catalogue to the IRAC-3.6 /im coordinates of the full SPIRE 
sample, within a 2 arcsec radius and find that 37 out of 84 
sources can be matched with a formal X-ray counterpart in 
at least one of the soft or hard X-ray bands (44 per cent). 

In order to reach higher completeness with respect to 
the X-ray association rate, we opt to go deeper by using 
the IRAC-3.6 )im position of each SPIRE source to extract 
counts from the X-ray map down to a Poisson significance 
< 3 x 10 -3 (corresponds to > 3<r in the case of a normal dis- 
tribution). For each object, the X-ray counts within the 70 
per cent encircled energy fraction (EEF) radius are summed 
up and a local value for the background is estimated using 
an annulus centered on the source with an inner radius 1.5 
times the 90 per cent EEF and a width of 50 arcsec (100 
pixels). For the estimation of the background, X-ray sources 
with Poisson significance < 4 x 10~ 6 are removed by exclud- 
ing pixels within the 95 per cent EEF radius of each source. 
The probability that the observed counts are a random fluc- 
tuation of the background is estimated using Poisson statis- 
tics. As in the IC catalogue, fluxes are estimated using a 
Bayesian methodology similar to that described in Laird et 
al. (2009) by integrating the source counts within the 90 
per cent EEF aperture. The count rates of sources in the 
soft and hard bands are converted to fluxes (ergs -1 cm~ 2 ) 
in the 0.5-2 and 2-10 keV energy intervals respectively by 
adopting a power-law X-ray spectrum with index T = 1.9, 
absorbed by the Galactic gas column density in the direction 
of the CDFN field, A?h=2x10 20 cm -2 . The uncertainties on 
the fluxes correspond to the 68 per cent confidence level. 
We choose T = 1.9 rather than the T — 1.4 used in the IC 
catalogue, as the former is consistent with the typically soft 
X-ray spectra of star-forming galaxies and unabsorbed AGN 
(e.g. Nandra & Pounds 1994) and hence more appropriate 
for our study. The X-ray detection rate now increases to 55 
per cent (46 out of 84 objects have an X-ray flux of > 3<r in 
at least one X-ray band). For objects without a formal or low 
significance X-ray counterpart, the 99.7 per cent confidence 
upper limit for their fluxes is estimated. 

We also estimate hardness ratios (HRs), defined as (H — 
S)/(H + S) where H and S refer to the hard (2-7keV) 
and soft (0.5-2 keV) count rates (counts s -1 ). The hardness 
ratio describes the relative strength of the hard to soft band 
signal and hence gives an indication of the shape of the X-ray 
spectrum. The X-ray properties of the sample are displayed 
in tabled 



2.3 The SPIRE sample 

In the full SPIRE sample, 67 out of 84 sources (80 per cent) 
are associated with a spectroscopic redshift and 16 with 
a photometric redshift only (1 has neither a spectroscopic 
nor a photometric redshift) — see Raymond et al. in prep.; 
also Eales et al. 2010 for photometric redshifts. The spec- 
troscopic redshifts come from the Barger, Cowie & Wang 
(2008) spectroscopic survey of the GOODS-N 24 fjxn popu- 
lation, the survey being more than 90 per cent complete at 
f24/jm >250/^Jy. We investigate the reliability of the photo- 
metric redshifts and find that although there is good agree- 
ment below z~l, most of the z >1 photometric redshifts 
significantly deviate (Az >0.5) from the true redshift (see 
Fig. [l}. Hence, by extrapolation, the redshift estimates for 
the majority of our photometric sub-sample (11 objects with 
£phot >1) are unlikely to be reliable. We therefore choose to 
restrict the work presented in this paper to the 67 spectro- 
scopically identified SPIRE sources. This excludes 4 X-ray 
detected sources, only 2 of which are potentially relevant to 
our study, the others having hard hardness ratios character- 
istic of AGN (see section |2T4"1 on the selection of AGN). As a 
result, the 67 sources with spectroscopic redshifts constitute 
a near-complete sample for our purposes and hereafter, for 
simplicity, we refer to them as as the 'SPIRE sample'. 

The SPIRE sample of 67 sources spans a redshift range 
of 0.2 < z < 2.6, with a mean of (z) =0.86. We calculate total 
infrared luminosities (Lir/L©, 8-1000 (an), by performing a 
X 2 fit on the available photometry (24-500 ^m), using the 
entire range of the Siebenmorgen & Kriigel (2007, SK07) li- 
brary which consists of about 7000 templates. In Symeonidis 
et al. (2008) we showed that the SK07 library can accommo- 
date a large range of SEDs and multiple dust temperature 
components and is thus ideal for the dusty starburst sam- 
ple we have selected here. Following the method outlined 
in Symeonidis et al. (2009), we calculate the total infrared 
luminosity (Lib,, 8-1000 /jm) and corresponding uncertainty 
for each object. The resulting lcr uncertainty in Lir (typi- 
cally of the order of 5 per cent) is combined with a systematic 
uncertainty of 8 per cent which takes into account the aver- 
age calibration uncertainty of MIPS and SPIRE. The final 
fractional uncertainties on Lib. are comparable to the frac- 
tional uncertainties on the SPIRE fluxes. We find 38 sources 
to be in the luminous infrared galaxy regime (LIRGs, 57 per 
cent), 22 are ultraluminous infrared galaxies (ULIRGs, 33 
per cent) and 7 are normal IR galaxies (Zir/Lq < 10 11 ). 

To calculate X-ray luminosities, the soft and hard fluxes 
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Figure 1. Photometric versus spectroscopic redshifts for the full 
SPIRE sample of 84 sources. The solid line represents a one-to- 
one relationship. The 43 X-ray detected objects are denoted with 
green diamonds. All objects at (x,y) = (0,p) have no spectro- 
scopic redshift and a photometric redshift of p. Note that most 
photometric redshifts above z~l, deviate significantly from the 
true redshift. 



are ii-corrected using a photon index F=1.9; hereafter Lsx 
and L H x refer to the soft (0.5-2 keV) and hard (2-10 keV) 
X-ray luminosities. For objects with X-ray upper limits, we 
use stacking analysis (described in Georgakakis et al. 2008) 
to explore their mean X-ray flux, by grouping them into 
LIRGs and ULIRGs and stacking separately in the soft and 
hard X-ray bands; see table [T] for the stacking results. For 
the normal IR galaxies we are unable to perform stacking as 
only 2 sources are not X-ray detected. X-ray luminosities for 
the stacked groups are computed using the mean redshift of 
the group. 

Fig. [2] shows the X-ray/IR luminosities of the SPIRE 
sample, with the AGN-classified sources clearly marked on 
the plot (see section 12.41 for details on AGN classification) . 
Note that due to the lower sensitivity of the hard X-ray 
band, many sources have no hard X-ray counterpart — the 
soft X-ray association rate is 54 per cent whereas the hard 
X-ray association rate is 29 per cent. 



2.4 The AGN fraction 

As our analysis focuses on the X-ray/IR properties of star- 
forming galaxies, we require the candidate sources to be 
dominated by star-formation in both the X-rays and the in- 
frared. The latter is ensured as our sample is submm-bright; 
as mentioned earlier, the dust masses implied for galaxies 
whose SEDs peak in the far-IR/submm, as well as their high 
energy output in that part of the spectrum, suggests that 
infrared emission is strongly linked to star-formation. Weed- 
ing out AGN from the sample therefore rests on diagnosing 
the origin of the X-ray emission. 

We investigate the incidence of AGN with the aid of 2 




Figure 2. The X-ray versus IR luminosities of the 67 SPIRE 
sources with redshifts — on the top panel we plot Lgx vs I/tr 
and on the lower panel Lhx vs Ltr. X-ray detected sources are 
denoted by black filled circles. Non-detections are incorporated 
into the stacked groups denoted by the 2 large stars (see table 
[T] for details on the stacking). Note that for the ULIRG cate- 
gory there is no stacked signal in the hard band and we plot the 
3 it upper limit. AGN classified objects are marked on the plot 
by green triangles, red diamonds and blue squares, depending 
on the criterion used to identify them. Green triangles are hard 
sources with HR>-0.1, blue squares correspond to AGN identi- 
fied through high excitation or broad lines in their optical spectra 
and red diamonds represent log [fHx/frt] >-l sources. Note that 
all, bar one, log [fHx/fR.] > -1 sources also have a hard X-ray lumi- 
nosity greater than 10 42 ergs" 1 . The two LIRGs and one ULIRG 
which are hard X-ray detected, but have no associated AGN sym- 
bol, are ambiguous classifications (see Section [2.4| l. Although they 
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X-ray criteria: the hard X-ray to R-band flux ratio (rkx/fR,) 
and the hardness ratio, relying on the well established re- 
sult that AGN are the most powerful X-ray emitters and 
that hard X-rays are minimally absorbed by the torus gas 
and dust. The fx/frt ratio has been extensively used to sepa- 
rate AGN and starburst systems since early observations of 
spectroscopically identified AGN have shown them to domi- 
nate the -l<log [fx/fii]<l parameter space (e.g. Schmidt et 
al. 1998; Akiyama et al. 2000), with star- forming galaxies 
having values of log [fx/fi)]<-l and typically <-2 (e.g. Horn- 
schemeier et al. 2001, 2003; Akiyama et al. 2003; Georgan- 
topoulos, Georgakakis & Koulouridis 2005; Georgakakis et 
al. 2007) . Our second criterion (the hardness ratio) is related 
to the differences between observed X-ray spectra of AGN 
and star-forming galaxies. Although both star-forming sys- 
tems and unabsorbed AGN are characterised by soft X-ray 
spectra and hence low hardness ratios, absorbed AGN have 
hard spectra and high values of hardness ratio, as low X- 
ray frequencies are more severely attenuated by gas column 
densities in the line of sight. Measured photon indices (r) of 
starburst galaxies range between 3 and 1.2 (e.g. Franceschini 
et al. 2003; Lehmer et al. 2010), corresponding to hardness 
ratios between -0.65 and -0.1 and therefore a cut at HR=-0.1 
should successfully isolate absorbed AGN. 

Although the properties of star-forming galaxies have 
not been seen to deviate from log [fx/frt];$-l, low luminosity 
AGN do exhibit properties which overlap with the starburst 
parameter space and hence will be missed. In addition, as 
this indicator is linked to a hard band detection, AGN which 
are only soft-band detected will not be identified. Finally, 
the hardness ratio diagnostic neglects all unabsorbed AGN. 
Aiming for a higher completeness in the AGN fraction, we 
supplement our diagnosis with optical spectroscopy, through 
inspection of the spectral classification of the X-ray-detected 
SPIRE sources in Trouille et al. (2008). In total, we find 25 
X-ray sources with an assigned spectral class as defined in 
Barger et al. (2005); 5 are high-excitation-line AGN (sources 
with [NeV] or CIV or strong [OIII] lines, i.e. EW([OIII] 
A 5007) > 3EW(H/3), 1 is a broad line AGN (FWHM line 
widths greater than 2000 km s -1 ) and the rest are classed 
as star-forming galaxies (strong Balmer lines and no broad 
or high-ionization lines) — see table [2] for a summary. An 
X-ray-detected SPIRE source which satisfies at least one 
of: log [fHx/fii]>-l, HR>-0.1 and high-excitation or broad 
spectral lines, is classed as AGN-dominated in the X-rays 
and hence not included in our study of the X-ray/IR star- 
formation correlation. As we shall see below, the combina- 
tion of these diagnostics enables us to pinpoint all powerful 
AGN in the sample, as well as identify some less luminous 
ones. 

Fig. [2] shows that most hard X-ray detected objects are 
classed as AGN, whereas due to the higher sensitivity of the 
soft-band, the soft X-ray detections are an almost even mix- 
ture of star-forming galaxies (SFGs) and AGN. In addition, 
as stated above, it is likely that some of SFGs host AGN 
not identified by our criteria. There are also 3 sources (ob- 
jects 35, 43 and 60; see caption of Fig. [2] and table [2J) which 
are only marginally classified as AGN, because of highly 
uncertain or unconstrained HR and/or fHx/fit values. We 
incorporate them in the star-forming sub-sample, although 
our conclusions do not change if we instead include them 
in the AGN sub-sample. It is also worth pointing out, that 



all, bar one, log [fHx/frt]>-l sources have a hard-band X-ray 
luminosity greater than 10 42 ergs _1 . An X-ray luminosity 
cut at 10 42 ergs _1 is widely used in extragalactic surveys 
to separate AGN and star-forming galaxies (e.g. see reviews 
by Fabbiano 1989; Brandt & Hasinger 2005), based on the 
fact that in the local Universe, where high resolution X-ray 
studies are possible, the most X-ray luminous starburst sys- 
tems have X-ray luminosities up to 10 42 ergs _1 in the full 
band and up to a few times 10 41 ergs -1 in the hard band 
(see Moran, Halpern & Helfand 1996, 1999; Zezas, Georgan- 
topoulos & Ward 1998). Although, we do not explicitly use a 
luminosity cut here, so as not to bias our sample against X- 
ray luminous star-forming sources, we see that its use would 
produce entirely consistent results with the AGN identified 
through the fex/fR ratio. 

In total we class 18 objects as AGN-dominated in the 
X-rays, a fraction of 27 (±10) per cent. This is higher than 
the X-ray AGN fraction of the Spitzer 70 /im-selected sam- 
ples of Symeonidis et al. (2010) and Kartaltepe et al. (2010) 
who find only 7 and 9 per cent X-ray-selected AGN respec- 
tively. We attribute these differences to (i) the much shal- 
lower X-ray surveys that their results are based on and (ii) 
differences in the relative fractions of LIRGs and ULIRGs 
in the 3 samples. With respect to the first point, the 2Ms 
Chandra survey that we use here implies a large increase in 
on-axis sensitivity from a 250 ks survey, so some of the lower 
luminosity AGN which are picked up with the 2 Ms survey 
would be missed by a shallower X-ray survey. With respect 
to the second point, our SPIRE sample contains a higher 
fraction of ULIRGs, which would imply an overall higher 
AGN fraction, as various observations point towards an in- 
crease of the AGN fraction as a function of Lir (e.g. Veilleux 
et al. 1995; 1999; Lutz et al. 1998; see also Symeonidis et al. 
2010; Kartalepe et al. 2010). 

On the other hand, the X-ray AGN fraction of the 
SPIRE sources is lower than the Alexander et al. (2005b) 
work on the X-ray properties of SMGs. They find about 
75 per cent of their sample hosting X-ray AGN, although 
as they point out, there is an additional bias due to the 
radio (1.4 GHz) selection, resulting in the best estimate of 
the AGN fraction being closer to 38 per cent (Alexander et 
al. 2005a). Recent studies by Laird et al. (2010) and Geor- 
gantopoulos, Rovilos & Comastri (2011) of purely submm- 
selected SMGs show X-ray AGN fractions of 20-29±7 per 
cent and 29±9 per cent respectively, entirely consistent with 
our SPIRE AGN fraction. 



3 THE SOFT X-RAY/IR STAR-FORMATION 
RELATION 

3.1 Trends in the Lsx/^m ratio 

In order to examine the X-ray/IR properties of the SPIRE 
SFGs, we examine the Lsx/im ratio as a function of Lir 
in Fig. [3] For completeness, we also show the sources host- 
ing AGN, clearly marked on the plot. However we do not 
include them in our analysis of the X-ray/IR relation, focus- 
ing solely on star-forming galaxies. Note that, due to small 
number statistics, any kind of evolution in the Lsx/Lir ratio 
cannot be measured within the SPIRE SFGs. Nevertheless, 
we opt to compare their properties with other samples taken 
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from the literature. These are the local IRAS 60 ^m-selected 
galaxies from David, Jones & Forman (1992), the local LIRG 
sample from Lehmer et al. (2010; soft X-ray data are cour- 
tesy of B. Lehmer), the local LIRG sample from Iwasawa et 
al. (2009; 2011), local/low-redshift ULIRGs from Teng et al. 
(2005), local ULIRGs from Franceschini et al. (2003), local 
ULIRGs from Grimes et al. (2005), local optically identified 
HII galaxies and low redshift X-ray/radio identified star- 
forming galaxies from Ranalli, Comastri & Setti (2003), low 
redshift X-ray selected CDFS sources from Rosa-Gonzalez et 
al. (2007) and moderate redshift ((z)=0.6), 70 ^m-selected 
galaxies from Symeonidis et al. (2010). In order to be con- 
sistent, we base our comparison on objects from these sam- 
ples, which do not host identifiable AGN (AGN are iden- 
tified by various methods in these studies) and with a far- 
IR (40-120 (im) or total IR (8-1000 /im) luminosity in the 
10< log (Lir/L©) <13 range. Where required, we make ap- 
propriate conversions to the cosmology we use in this paper, 
as well as to the soft X-ray band range of 0.5-2 keV. Wher- 
ever necessary we also convert far-IR luminosity (Lfir) to 
total infrared luminosity (Lir) assuming Lfir is 70 per cent 
of Lib. (Helou et al. 1988). All X-ray luminosities are in the 
0.5-2 keV band, corrected only for Galactic absorption and 
L"-corrected (with T=1.9). 

In Fig. [3] we also plot the linear relation from Ranalli, 
Comastri & Setti (2003), as above appropriately modified by 
converting far-IR luminosity (Lfir) to total infrared lumi- 
nosity (Lir). Note that the Ranalli et al. relation was derived 
predominantly with a Lir <10 11 L© sample, a large fraction 
of which were normal star- forming galaxies (Lir <10 Lq) 
and is not necessarily representative of higher luminos- 
ity sources. However for comparison pusposes we extrap- 
olate it up to Lir = 10* Lq. Although Ranalli et al. discuss 
that their sample is consistent with both Lsx oc Lir and 
LsxocLjr 7 , our subsequent comparison uses the linear re- 
lationship as a point of reference, as it describes a scenario 
where X-ray emission scales with the star-formation rate. 
The linear Ranalli et al. relation derived for local galaxies 
is also seen to hold at low (z <0.4) redshifts for sources of 
equivalent luminosity, e.g. see Rosa-Gonzalez et al. (2007). 

Apart from the Symeonidis et al. (2010) sources, all 
other comparison samples are at z <0.4, with the majority 
at z <0.2. In order to quantify any change in the Lsx/Lir 
ratio as a function of IR luminosity, we group them in 5 in- 
frared luminosity bins (10< log (Lir/L©) <12.5 in steps of 
0.5 dex), calculating the mean Lsx/Lir ratio and la en- 
velopes for each bin, shown as a smooth grey shaded region 
in Fig. [21 Although the scatter is large, the average Lsx/Lir 
ratio of these 9 z <0.4 samples is seen to decrease with in- 
creasing IR luminosity: for log (Lir/Lq) >11.5 sources, the 
mean ratio is 0.7 dex lower than for normal IR galaxies. 
The linear Ranalli et al. relation (red line in Fig. |3J, which 
predicts a constant Lsx/Lir ratio of -3.84, overpredicts the 
properties of the more LiR-luminous sources. The Lsx/Lir 
ratio for the latter drops to an average value of -4.55±0.05 
(la scatter of 0.37). Taking the soft X-ray emission to be 
proportional to Lir (and hence the star- formation rate), 
the Lsx/Lir relation for local Lir >3 xlO 11 L0 galaxies can 
thus be written as: 

log L sx = log Lir - 4.55 (1) 

The result that more LiR-luminous sources are soft X- 



ray deficient with respect to their IR emission compared to 
their less luminous counterparts has been noted before by 
Franceschini et al. (2003) and Grimes et al. (2005). Grimes 
et al. (2005) discuss that although the X-ray spectra of both 
normal IR galaxies and ULIRGs are dominated by the ther- 
mal component at low energies < 2 keV, the latter have a 
much larger characteristic size of their diffuse X-ray emit- 
ting regions; the radius which encloses 90 per cent of the flux 
extends up to 30kpc in ULIRGs. However, the discrepancies 
in soft-band X-ray emission cannot be reconciled by missing 
flux, as Grimes et al. show that, once bright point source flux 
is removed, the surface brightness at the 90 per cent radius 
is roughly constant in both galaxy types. Instead, the obser- 
vation that the diffuse X-ray emission in ULIRGs extends to 
much larger scales suggests that the density and hence emis- 
sivity of the X-ray emitting gas might be lower, which could 
in part be responsible for making high luminosity sources 
soft-X-ray-deficient . 

The low redshift samples we use here are highly com- 
plete, with the X-ray association rate greater than 80 per 
cent, higher than the association rate of 51 per cent for 
the SPIRE SFGs (out of the 49 sources in our sample clas- 
sified as SFGs, only 25 are X-ray detected). This implies 
that for the comparison between these low redshift sam- 
ples and the SPIRE SFGs, it is appropriate to include the 
stacked results for the latter. Fig. [3] shows that the SPIRE 
SFGs have ratios consistent with the scatter seen in their 
local equivalents, within 2a. The average Lsx/Lir ratio for 
LIRGs and ULIRGs in the SPIRE sample is estimated using 
the detections and average stacked luminosities, the latter 
weighted by the number of sources in the stacked groups. We 
find (Lsx/Lir)= -4.43 with a la scatter of 0.37, consistent 
with the average ratio of their local counterparts (equation 
[TJ. Similarly, the normal IR galaxies (Lir<10 1 Lq) in the 
SPIRE sample appear to have properties consistent with 
those of local normal IR galaxies, described by the Ranalli 
et al. relation, although due to small number statistics this 
result is uncertain. The SPIRE AGN display large scatter 
in log [Lsx/Lir] and many clearly stand out with values 
>-4. However, some show ratios consistent with the SPIRE 
SFGs, either because their soft X-ray emission is intrinsically 
low or because of substantial absorption in the soft X-rays. 
Equivalently, some of the SFG-classified galaxies could in 
fact be AGN-dominated in the X-rays. Inevitably, separat- 
ing the AGN and SB contribution in these sources may only 
be possible with high resolution X-ray spectra, outside the 
capabilities of current X-ray observatories. 

Although the comparison between local and high red- 
shift SFGs is fairly straight forward, an issue which needs 
to be addressed is that at z>0, the soft band starts receiv- 
ing emission from higher energy X-rays. After ~2keV (rest- 
frame), the contribution from diffuse gas declines whereas 
that of point sources increases, with the X-ray spectrum 
changing from predominantly thermal to predominantly 
power-law (e.g. Persic & Rephaeli 2002; Franceschini et al. 
2003; Grimes et al. 2005). For most of the SPIRE SFGs this 
is not a serious problem since at z <1 the X-rays received by 
the soft band are still largely of thermal origin. On the other 
hand, for the higher redshift ULIRGs, the rest-frame 0.5- 
2keV part of the X-ray spectrum is mostly unconstrained 
by our data. The L"-correction with T=1.9 would be suffi- 
cient to correct for this, under the assumption that the X- 
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Figure 3. Lsx/^IR against Lm for the SPIRE SFGs (black filled circles) and AGN (black filled circles with red outline). The large filled 
stars correspond to the groups of stacked sources (sec table[l]for more details). For comparison, we also include various samples from the 
literature. Pink asterisks: the local IRAS 60 ^in-selected galaxies from David, Jones & Forman (1992), dark blue open circles: the local 
LIRG sample of Lehmer et al. (2010; soft X-ray data are courtesy of B. Lehmer), brown open squares: the local LIRG sample of Iwasawa 
et al. (2009; 2011), light blue open circles: local/low-redshift ULIRGs from Teng et al. (2005), green open diamonds: local ULIRGs from 
Franccschini et al. (2003), light green crosses: local ULIRGs from Grimes et al. (2005), small red open squares: local optically identified 
HII galaxies from Ranalli, Comastri & Setti (2003), large red open squares: low redshift X-ray/radio identified star-forming galaxies from 
Ranalli, Comastri & Setti (2003), dark blue plus signs: low redshift X-ray selected CDFS sources from Rosa-Gonzalez et al. (2007) and 
orange triangles (small for detections and large for stacking): moderate redshift ((z}=0.6), 70 /^m-selected galaxies from Symeonidis et 
al. (2010). The grey shaded region correponds to the smoothed mean and la envelopes of the local and low redshift samples calculated 
in 5 luminosity bins (10< log (Lir/Lq) <12.5 in steps of 0.5 dex). The vertical dotted lines separate the 3 IR luminosity classes that 
the sources fall into. The black horizontal line is equation [T] (solid part), extrapolated to lower luminosities than were used to derive it 
(dashed part). The red horizontal line corresponds to the Ranalli et al. (2003) linear soft X-ray/IR correlation (solid part; see Section 
I3.1J I extrapolated to higher luminosities than were used to derive it (dashed part). 
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ray SEDs of high redshift SFGs are similar to those of their 
local counterparts and hence relatively flat (in vi v ). This 
assumption, i.e. that the ratio of thermal to non-thermal 
emission stays relatively constant with redshift, is reason- 
able for two reasons: (i) z <1 SPIRE SFGs have Lsx/Lm 
ratios consistent with their local counterparts, so it would 
take fundamental changes in the origin of the soft X-ray 
emission for the z >1 ULIRGs in our sample to be different 
to z <1 LIRGs and ULIRGs; (ii) in the hard X-rays, where 
the 5-10 keV band receives emission from the non-thermal 
part of the spectrum at all redshifts, a simple if-correction is 
sufficient to rest-frame the flux and we find consistent prop- 
erties between local and high redshift sources (see section 
@J. To summarise: given that our data shows no evolution 
in the Lsx/Lm ratio up to z <1 and no evolution in the 
Lnx/im ratio up to z <2, it is reasonable to assume that 
the Lsx/Lm ratio is also constant up to z <2, although 
this would only be formally settled with the availability of 
<0.5 keV X-ray data. 



3.2 The Lsx/Lm ratio as a tracer of stellar 
evolution 

Mas-Hesse, Oti-Floranes & Cervifio (2008, hereafter 
MOC08) use the Lsx/Lm ratio as a tracer of stellar evo- 
lution and compute its variation over 30 Myr, based on 
the evolutionary synthesis models of Cervifio, Mas-Hesse & 
Kunth (2002), which assume solar metallicity and a Salpeter 
IMF with stellar mass limits of 2-120 M . In Fig. 0] we 
plot the MOC08 models (top panel) and the distribution 
in Lsx/Lm for the LIRGs and ULIRGs in the samples de- 
scribed in section [3TT1 (lower panel). According to MOC08, 
the Lsx/Lm ratio is highly sensitive to the star-formation 
history and the efficiency by which energy is coupled to the 
ISM and as a result it can span an order of magnitude at any 
given epoch. The MOC08 models shown here were derived 
for two different star-formation modes: stars forming dur- 
ing an instanteneous burst (IB) or an extend period of time 
(EB), although realistically, EB systems are not thought to 
be younger than ~10 Myr, as they are supposed to be a com- 
posite of many instantaneous bursts mimicing constant star 
formation. The MOC08 models are also shown as functions 
of e (=0.01 or 0.1), which is the fraction of energy incident 
onto the ISM that heats the X-ray emitting gas. Note that 
MOC08 neglect point sources, as these are estimated to con- 
tribute a negligible amount (up to 15 per cent) to the soft 
X-ray emission in the early stages of the starburst evolution 
examined by their models. 

MOC08 discuss that after a few Myrs, IR emission de- 
creases sharply in the IB mode as the most massive stars 
reach the end of their lifetimes and therefore the intensity 
by which the dust is being heated significantly decreases. On 
the other hand in the EB mode, IR emission stays relatively 
constant as massive stars are being continuously replenished. 
The soft X-ray luminosity undergoes a rapid increase during 
the first few Myrs, due to the injection of energy into the 
ISM from the winds of the most massive stars in both IB 
and EB modes. Cervifio, Mas-Hesse & Kunth (2002) pro- 
pose that, as a result, X-ray emission is strongly metallicity 
dependent during the first few Myr. Subsequently, injection 
of energy is dominated by supernovae so the soft X-ray lu- 
minosity increases for EB models, whereas it levels off in IB 
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Figure 4. In the top panel, we plot the Mas-Hesse, Oti-Florancs 
& Cervifio (2008) models which use the Lsx/Lm ratio as a tracer 
of stellar evolution and compute its variation over 30 Myr, based 
on the evolutionary synthesis models of Cervifio, Mas-Hesse & 
Kunth (2002), assuming solar metallicity and a Salpeter IMF with 
stellar mass limits of 2— 120 Mq. EB stands for extended bursts, 
IB stands for instantaneous burst, e is the efficiency by which me- 
chanical energy is used to heat up the X-ray emitting gas. In the 
lower panel we show the distribution in log [Lsx/Lm] ratio for the 
SPIRE LIRGs and ULIRGs (dashed and solid black histograms 
respectively) compared to equivalent sources from the samples de- 
scribed in section 13.11 These are local/low redshift LIRGs (pink 
region) and ULIRGs (green filled region) and the moderate red- 
shift log (Lir/Lq) >11 sources from Symeonidis et al. (2010; blue 
region). Note that for clarity, the comparison samples are not 
plotted as histograms, but rather we interpolate between the to- 
tal number of sources in each log [£sx/^ir] bin and shade the 
area underneath. Also, in the case of stacked data, each source 
in the stacked group is counted individually in the relevant bin. 
The log [Lsx/^ir] bins are 0.44 dex in width. 



models. When translated to a Lsx/Lm ratio, this implies a 
continuous increase in the ratio for IB models and a smaller 
rate of increase in the EB scenario. In addition, for a given 
efficiency and at any given epoch, the IB mode of star for- 
mation produces an X-ray excess compared to the EB mode, 
because of the larger number of supernovae that heat up the 
gas. 

In the lower panel of Fig. [4] we show the distribution 
in log [Lsx/Lm] for the SPIRE LIRGs and ULIRGs com- 
pared to the samples presented in section [3~T1 grouped into 
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local/low redshift LIRGs, local/low redshift ULIRGs and 
moderate redshift log(Lia/L0) >11 sources from Symeoni- 
dis et al. (2010). The log [Lsx/Lm] bins are 0.44 dex in width 
and in the case of stacked data, each source in the stacked 
group is counted individually in the relevant luminosity bin. 



As described in section 13.11 and also clearly seen here, 
there is significant overlap between the 5 groups, indicative 
of little or no change in the Lsx/Lir ratio of LIRGs and 
ULIRGs with redshift. However, note that stellar evolution 
might not be the only relevant factor in shaping the X-ray 
properties of SFGs. As discussed in Section 13.11 the char- 
acteristic size of the X-emitting region could influence the 
strength of X-ray emission and might explain why there is 
a group of local LIRGs with higher Lsx/Lir ratios than 
the bulk of the other objects (far-right section of the pink 
shaded region in the lower panel of Fig. ©) — these are the 
lowest luminosity LIRGs in the sample (Lm ~ 10 11 L©; see 
also Fig.©. 



Fig. [4] illustrates that the X-ray luminosities of the 
majority of objects are sufficiently low for the average 
log [Lsx/Lir] to be <-4 and according to the MOC08 
models, consistent with very young ages (<10Myr) and 
low efficiency (e=0.01). For the few remaining objects 
with log [Lsx/Lm] >-4, age, efficiency and mode of star- 
formation are degenerate, but apart from the EB, e =0.01 
model, they are also compatible with young ages up to 
10 Myr. Note that during the very early stages (<5 Myr) 
of star- formation, the Lsx/Lir ratio is sensitive to the IMF 
(slope, mass limits) — see Oti'-Floranes & Mas-Hesse (2010). 
For example, changing the IMF mass limits from 2-120 Mg 
in the MOC08 models to 0.1-100 M (e.g. Kennicutt 1998), 
can reduce the Lsx/Lm ratio by up to 30 per cent, whereas 
assuming a top-heavy IMF (slope= -1) would yield higher 
Lsx/Lir values in the IB scenario. Nevertheless, despite any 
changes in the IMF and/or efficiency (which cannot be quan- 
tified currently), the X-ray/IR ratio of LIRGs and ULIRGs 
is overall consistent with very young stellar ages. This might 
seem stange at first, as LIRGs and ULIRGs have high stel- 
lar masses and therefore a substantial older population. The 
fact that these systems are represented by young ages in 
the MOC08 models suggests that both the soft X-ray and 
infrared emission must be dominated by the current star- 
formation episode, with minimal contribution from older 
stellar populations or regions of lower SFR. This is not sur- 
prising as the main factors in heating the X-ray emitting gas 
are supernovae and stellar winds, both of which are charac- 
teristic of massive stars, implying that massive stellar evo- 
lution is the key player in the soft X-ray output of a high 
SFR galaxy. Moreover, although total infrared luminosity 
traces the starburst episode in LIRGs and ULIRGs up to 
~100Myr (e.g. see review by Kennicutt 1998 and references 
therein), many models (e.g. Devereux & Young 1990; Mas- 
Hesse & Kunth 1991; Cervino & Mas-Hesse 1994; Devereux 
& Hameed 1997) support that infrared luminosity peaks at 
< 10 Myr, suggesting that the peaks of IR and X-ray emis- 
sion are roughly coincident for galaxies undergoing starburst 
episodes. 



4 THE HARD X-RAY/IR STAR-FORMATION 
CORRELATION 

Fig- ©shows hard X-ray luminosity (Lhx) versus Lir for the 
SPIRE sample. This is also plotted as a ratio in Fig. © to 
facilitate comparison with Fig. © We also plot the Ranalli 
et al. Lhx-Lir linear relation, again appropriately modified 
to convert Lfir to Lm and extrapolated to higher lumi- 
nosities than were used to derive it (see section©. Some 
of the samples shown in Fig. © are also included here — 
the Ranalli et al. (2003) star-forming galaxies, the Frances- 
chini et al. (2003) ULIRGs, the Symeonidis et al. (2010) 
sample (consisting of stacked data only in the hard X-rays), 
the local LIRG/ULIRG samples from Iwasawa et al. (2009) 
and Lehmer et al. (2010) and the 2 < 0.16 ULIRG sample 
from Teng et al. (2005). Also included is the subnim galaxy 
(SMG) sample of Laird et al. (2010). All comparison samples 
comprise sources which do not host AGN according the clas- 
sifications given by the authors, however note that as for the 
soft X-rays, there is always a non- negligible possibility that 
some have been misclassified. Finally, all necessary conver- 
sions were performed to the data in order for the comparison 
to be consistent (see Section ©for an account). 

Figs. © and © show that local/low redshift 
LIRGs/ULIRGs have significantly lower hard X-ray 
luminosities than what is predicted by the Ranalli et al. 
(2003) relation derived for normal star-forming galaxies, 
a result that has been well established by various studies 
e.g. Persic & Rephaeli (2002), Franceschini et al. (2003), 
Persic et al. (2004), Iwasawa et al. (2009), Lehmer et al. 
(2010). This is consistent with the picture seen in the soft 
X-rays (Fig. ©, where the discrepancies are thought to be 
primarily due to differences in the extent of X-ray emitting 
region, the efficiency of energy coupling to the ISM and 
mode of star-formation. With respect to the hard X-rays, 
Persic & Rephaeli (2002), Franceschini et al. (2003) and 
Persic et al. (2004) argue that the differences in X-ray 
emission could arise because of varying contributions from 
HMXBs and LMXBs. In LIRGs and ULIRGs, emission 
from HMXBs dominates the hard X-ray luminosity, but 
in lower luminosity star-forming galaxies, LMXBs also 
contribute substantially, implying that for such sources, 
hard X-rays are not a reliable tracer of star-formation. 
Lehmer et al. (2010), derive an expression for the hard 
X-ray luminosity as a function of both the star-formation 
rate (SFR, i.e. emission from HMXBs) and the stellar mass 
(M, i.e. emission LMXBs), 



L H x = aM + /3SFR. 



(2) 



They point out that in the low SFR regime 
(SFR<5Mq yr _1 ), there seems to be an intrisic corre- 
lation between SFR and stellar mass. This causes the first 
term in equation © to dominate, resulting in the observed 
linear relation between Lhx and SFR (e.g. in the Ranalli 
et al. 2003 sample), even though the former quantity is 
not directly related to the latter; this is also shown in 
Grimm, Gilfanov & Sunyaev (2003); Gilfanov, Grimm & 
Sunyaev (2004); Colbert et al. (2004). In the high SFR 
regime (SFR>5M0 yr _1 ), Lehmer et al. show that stellar 
mass only weakly correlates with the star-formation rate 
[log (M) = 10.4+0.3 log (SFR)], so the relationship between 
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Figure 5. Hard X-ray luminosity (Lhx) versus Ltr for the SPIRE SFGs (black filled circles) and AGN (black filled circles with red 
outline); large filled stars correspond to the stacking groups (see table [l}. Note that only LIRG group has a stacked signal in the hard 
X-rays, whereas for the ULIRGs we plot the 3cr upper limit. For comparison we show the Ranalli et al. -Lhx-^ir linear relation (red 
solid line) extrapolated to higher luminosities than were used to derive it (dashed red line). Some of the samples shown in Fig. [3] are also 
included — the Ranalli et al. (2003) star- forming galaxies (red empty squares), the Franccschini et al. (2003) ULIRGs (green diamonds), 
the Symconidis et al. (2010) sample (orange triangles, consisting of stacked data only in the hard X-rays), the local LIRG/ULIRG samples 
from Iwasawa et al. (2009, brown squares) and Lehmer et al. (2010, blue circles) and the z < 0.16 ULIRG sample from Teng et al. (2005, 
blue circles). Also included is the submm galaxy (SMG) sample of Laird et al. (2010; pink crosses). The solid black line is the Lehmer et 
al. (2010) equation (equation [3] in Section|4] equivalent to equation [4] for Ljr >10 Lq) and 2 a scatter of 0.68 dex (dashed black lines). 



Lhx and Ltr becomes 



L HX = (1.9 x 10 26 )L? H 3 + (4.15 x l(T 5 )Li R 



where SFR(M Q yr- 1 ) 
niccutt et al. 1998), 



/3=1.62xl0 39 ergs _1 (M 



©yr 



9.8 x KT^Lir (Ken- 
=9.05 xlO 28 ergs" 1 M" 1 and 
)~ 1 and both luminosi- 



ties are in ergs . This is plotted in Figs. [5] and [6] together 
with the 2<r scatter of 0.68 dex quoted in Lehmer et al. 
(2010). Note that for high SFRs and L m >10" L a simpler 
version of this equation is derived by dropping the first 
term in equation [2 whereby Lhx becomes proportional to 
Lip.: 



log Lhx = log Ltr - 4.38 



This is about 0.6 dex lower than the relation of Ranalli et 
al. (2003) derived for normal star-forming galaxies. 

(3) Note that due to the lower sensitivity of Chandra in 
the hard X-rays (almost an order of magnitude lower than 
in the soft X-rays), we only detect 4 SPIRE SFGs, albeit at 
low signal to noise, 3 of which are ambiguous and could po- 
tentially be AGN dominated (see Section ^. 4p . Nevertheless, 
our overall results, which include the stacking, are in agree- 
ment with the Lehmer et al. relation (equation [4]) derived 
for local LIRGs and ULIRGs within the 2<r scatter. This 
is also the case for the moderate and high redshift samples 
of Symeonidis et al. (2010) and Laird et al. (2010), on par 
with a no-evolution scenario for the Lhx/Lir ratio of SB- 
dominated LIRGs and ULIRGs up to z~3. In terms of the 

(4) SPIRE AGN, we see that most occupy a different parameter 
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Figure 6. The Lhx/^ir ratio plotted against Lir for the SPIRE 
SFGs (black filled circles) and AGN (black filled circles with red 
outline); large filled stars correspond to the stacking groups (see 
table [TJ- The vertical dotted lines separate the 3 IR luminosity 
classes that the sources fall into. All other symbols and lines are 
the same as in Fig. [5] 



high-excitation or broad optical lines. Out of the 43 X-ray 
detected sources, 18 are identified as AGN, amounting to a 
fraction of 27 (±10) per cent in the SPIRE sample. Note 
that our aim is to exclude any obvious AGN, but we cannot 
ensure a pure star-forming sample, since low luminosity 
and/or unabsorbed AGN will be missed. Once obvious 
AGN are eliminated, our analysis is based on 24 X-ray 
detected SPIRE objects and 23 non-X-ray detected SPIRE 
sources, whose X-ray properties we assess through stacking. 

We conclude that the X-ray /IR properties of star- 
formation-dominated SPIRE galaxies, quantified by the 
soft and hard X-ray/IR luminosity ratios (Lsx/Ltr, and 
Lhx/Lir), are consistent with those observed in the local 
(z < 0.1) Universe. They are also consistent with the ra- 
tios of equivalent sources at more moderate (z ~ 0.6) and 
higher redshifts [z ~ 2), indicating that, overall, there is no 
evidence for evolution of the X-ray/IR correlation for IR- 
luminous galaxies with redshift. Note that, although none 
of the samples we present here are entirely impervious to 
AGN contamination, the main outcome of our study is ro- 
bust, in the sense that once all identifiable bona-fide AGN 
are excluded, the remaining sources have consistent average 
X-ray/IR ratios at all redshifts. 

The linear nature of the X-ray/IR correlation describes 
a scenario whereby both X-ray and IR emission trace the 
star-formation rate and is thus important in calibrating X- 
ray emission as a star formation indicator, providing an al- 
ternative census of the global SFR density. This will be par- 
ticularly applicable in future deeper X-ray surveys, where 
the X-ray number counts at the faint flux end will be dom- 
inated by star-forming galaxies. 



space to the star-forming galaxies. This suggests that in the 
absense of hardness ratio, R-band magnitude or optical line 
measurements, the hard X-ray/IR correlation can success- 
fully be used to separate AGN from star-forming galaxies. 



5 SUMMARY AND CONCLUSIONS 



Our analysis is based on a sample of 67, 
spectroscopically-identified, SPIRE-detected, submm- 
bright (/250 M m >17.4mJy; >3<r including confusion noise) 
galaxies in the GOODS-N field, observed as part of the 
H ers chel /HerMES survey, the majority of which are LIRGs 
and ULIRGs at z ~1. Taking advantage of the 2 Ms Chandra 
Deep Field North survey, we investigate the SPIRE sources' 
soft and hard X-ray emission relative to their infrared 
output and compare our results to several other samples 
from the literature. The focus of our study is the X-ray/IR 
correlation for star- forming IR- luminous galaxies at z~l. 
Our sample is selected to be submm-bright, implying that 
all objects are star-formation-dominated in the infrared, 
thus we identify AGN through X-ray diagnostics and optical 
spectral lines. We class as AGN-dominated in the X-rays, 
objects which satisfy at least one of the following criteria: a 
hard X-ray/R-mag flux ratio (log [fax/fR.]) greater than -1, 
a hardness ratio (HR) greater than -0.1 and the presence of 
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Table 2. The full SPIRE sample of 84 SPIRE-detected, submm-bright (/ 2 50 M m >17.4mjy) galaxies in the GOODS-N field, observed as 
part of the Herschel/HerMES survey. The columns are as follows: (1) ID, (2) (3) IRAC-3.6^m RA and DEC, (4) spectroscopic redshift 
(available only for 67 sources), (5) total infrared luminosity (8-1000 /im), (6) the cross- matched ID from the IC catalogue of HDFN 
Chandra data, (7) (8) (9) (10) the flux and luminosity in the soft and hard bands (SX stands for 'soft X-ray band' and HX stands for 
'hard X-ray band'), (11) the hardness ratio (HR), (12) the /hx//r ratio used to identify AGN (see section [23J , computed for sources 
which are detected in both the R-band and the hard X-rays. Note that, although we provide the corresponding HDFN X-ray ID from 
the IC catalogue, the fluxes and upper limits we present here are all derived using the IRAC-3.6/xm position of each SPIRE source to 
extract counts from the X-ray map down to a Poisson significance < 3 X 10 -3 — see Section l2.2l for details. '< fx' corresponds to the 3cr 
upper limit. Column (13) refers to the classification of the source by optical spectroscopy (Barger et al. 2005; Trouille et al. 2008). SFG 
stands for 'Star Forming Galaxy', sources with strong Balmer lines and no broad or high-ionization lines. HEL AGN stands for 'High 
Excitation Line AGN', sources with [NeV] or CIV lines or strong [OIII] lines, i.e. EWQOIII] A5007) > 3EW(H/3). BL AGN stands for 
'Broad-line AGN', sources with FWHM line widths greater than 2000 km s -1 . Note that the Trouille et al. (2008) classification applies 
only to X-ray sources. Where there is no information available, the spectrum could not be classified. Column (14) refers to whether 
the source is SB or AGN dominated in the X-rays, according to the criteria outlined in section 12.41 The 3 ambiguous sources (labelled 
'Amb'), are included in the SFG sub-sample, but our conclusions do not change if instead they are included in the AGN sub-sample. 



ID 


RA 


Dec. 


z 


^IR 


HDFN ID 


/sx 


/hx 


^SX 


ifiX 


HR 


/hx//r 


Spectral 


X-ray 




deg 


deg 




l°g(L©) 




erg/s/cm 2 


erg/s/cm 2 


log (erg/s) 


log (erg/s) 






Class 


Class 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


1 


189.28772 


62.38866 


1.26 


12.39 


— 


<2.37e-16 


<7.67e-16 


— 


— 


— 


— 


— 


— 


2 


189.27887 


62.38520 


0.50 


11.42 


— 


<1.39e-16 


<2.98e-16 


— 


— 


— 


— 


— 


— 


3 


189.29427 


62.37627 


1.52 


12.39 


— 


<1.34e-16 


<3.39e-16 


— 


— 


— 


— 


— 


— 


4 


189.27174 


62.33079 


0.84 


11.75 


— 


<7.65e-17 


<8.15e-16 


— 


— 


— 


— 


— 


— 


5 


189.35646 


62.32817 


0.28 


10.75 


hdfn-358 


1.86e-16 


<3.93e-16 


40.64 


— 


-1 


— 


SFG 


SFG 


6 


189.14424 


62.32388 


0.87 


11.72 


hdfn-086 


1.62e-16 


9.04e-15 


41.75 


43.50 


83 - 86 


60+ 003 


SFG 


AGN 


7 


189.08646 


62.31671 


— 


— 


— 


<5.26e-17 


<2.79c-16 


— 


— 


— 


— 


— 


— 


8 


189.54642 


62.30510 


— 


— 


— 


<2.44e-16 


<7.17e-16 


— 


— 


— 





— 


— 


9 


189.22009 


62.30213 


0.25 


10.90 


— 


2.05e-17 


<1.94e-16 


39.59 


— 


-1 





— 


SFG 


10 


189.44516 


62.29540 


1.15 


11.81 


— 


<l.llc-16 


<3.15e-16 


— 


— 


— 


— 


— 


— 


11 


189.39382 


62.28980 


0.64 


11.62 


— 


4.24c-17 


<1.24e-16 


40.85 


— 


-1 





— 


SFG 


12 


189.51512 


62.28651 


0.28 


11.01 


— 


<1.55e-16 


<3.56e-16 


— 


— 


— 





— 


— 


13 


189.08741 


62.28600 


1.99 


12.63 


— 


<7.07e-17 


<3.02c-16 


— 


— 


— 





— 


— 


11 


189.57885 


62.28439 


— 


— 


— 


<7.20e-17 


<6.14e-16 


— 


— 


— 





— 


— 


15 


189.13544 


62.28315 


0.44 


11.26 


— 


<2.61e-17 


<2.75e-16 


— 


— 


— 


— 


— 


— 


16 


189.31910 


62.27872 


0.56 


11.28 


hdfn-264 


4.54c-17 


1.95e-16 


40.73 


41.37 


U.IM_q 20 


-, ,0+0.22 
-l-tS-0 .78 


HEL AGN 


AGN 


17 


189.37984 


62.27229 


0.98 


11.89 


— 


<7.81e-17 


<1.31c-16 


— 


— 


— 





— 


— 


18 


188.99886 


62.26384 


0.38 


11.38 


hdfn-335 


1.55e-16 


<4.42e-16 


40.86 


— 


-1 





SFG 


SFG 


19 


189.07638 


62.26402 


2.00 


12.70 


hdfn-291 


3.01c-17 


<3.87c-16 


41.89 


— 


-1 





— 


SFG 


20 


188.99144 


62.26019 


0.38 


11.60 


hdfn-503 


4.66e-17 


<4.47e-16 


40.34 


— 


-1 





— 


SFG 


21 


189.06708 


62.25379 


2.58 


12.91 


hdfn-117 


2.01c-16 


1.12e-15 


42.98 


43.72 


24+ 008 


69+ 009 

u - D9 -o.io 


— 


AGN 


22 


189.48220 


62.25186 


0.19 


10.53 


hdfn-356 


4.61c-16 


5.48e-18 


40.66 


38.74 


W+ 013 

-U.OM_ 14 


q opr + 1.26 

-o-OO_0.56 


— 


SFG 


23 


188.96211 


62.24923 


0.39 


11.03 


— 


<2.65e-16 


<1.41e-15 


— 


— 





— 


— 


— 


21 


189.14271 


62.24249 


0.52 


11.20 


— 


<2.49e-17 


<1.77c-16 


— 


— 





— 


— 


— 


25 


189.14832 


62.23999 


2.01 


12.88 


hdfn-082 


1.56c-16 


1.73e-15 


42.61 


43.65 


u -°' -0.07 


25+ ' 05 
u - zo -0.05 


HEL AGN 


AGN 


26 


189.53851 


62.23786 


— 


— 


— 


<1.08e-16 


<7.02e-16 


— 


— 


— 





— 


— 


27 


189.26141 


62.23376 


1.25 


12.03 


— 


<3.21e-17 


<1.14e-16 


— 


— 


— 


— 


— 


— 


28 


189.15212 


62.22800 


0.56 


11.09 


— 


<2.06e-17 


<1.08e-16 


— 


— 





— 


— 


— 


29 


188.97445 


62.22703 


0.88 


12.00 


— 


<6.34e-17 


<3.07c-16 


— 


— 


— 


— 


— 


— 


30 


189.50379 


62.22669 


0.44 


11.77 


hdfn-003 


3.72e-14 


4.19e-14 


43.41 


43.46 


-0 60+ 001 
u - DU -o.oi 


-0 31+ 001 


BL AGN 


AGN 


31 


189.29727 


62.22530 


2.00 


12.82 


hdfn-173 


5.75c-17 


3.56e-16 


42.17 


42.96 


11+ al5 
U - ia -0.14 


_ n4 c:+0.11 


— 


AGN 


32 


189.37812 


62.21623 


— 


— 


— 


<4.86e-17 


<1.43e-16 


— 


— 






— 


— 


33 


189.08131 


62.21457 


0.47 


11.52 


hdfn-112 


8.15c-17 


<5.19e-16 


40.82 


— 


-1 


— 


HEL AGN 


AGN 


31 


189.42242 


62.21422 


1.60 


12.52 


hdfn-408 


5.57c-17 


5.03e-16 


41.93 


42.88 


o 7 +o.i9 

u -°'-0.18 


-o.3oi° : g 


SFG 


AGN 


35 


189.14381 


62.21139 


1.22 


12.56 


hdfn-087 


6.15c-17 


1.87e-16 


41.69 


42.17 


-o.23± ;g 


-0.861°;- 


SFG 


Amb 


36 


189.14571 


62.20672 


0.56 


11.30 


— 


<3.01e-17 


<9.47e-17 


— 


— 






— 


— 


37 


189.42151 


62.20581 


— 


— 


— 


<5.30e-17 


<1.83e-16 


— 


— 


— 


— 


— 


— 


38 


189.43598 


62.20522 


0.23 


10.41 


— 


<1.56e-16 


<1.65e-16 


— 


— 


— 


— 


— 


— 


39 


189.14363 


62.20360 


0.46 


11.65 


hdfn-088 


1.47e-16 


2.53e-16 


41.04 


41.28 


-0.38±8;i2 


2 10+ 012 

-z.ru_ 14 


HEL AGN 


AGN 


10 


189.40440 


62.20136 


0.41 


11.10 


hdfn-493 


8.42e-17 


<3.17e-16 


40.69 


— 


-1 


— 


— 


SFG 


11 


189.42067 


62.20013 


1.17 


12.06 


hdfn-008 


7.26e-15 


1.40e-14 


43.71 


44.00 


_n 39+0-02 


71+ 002 
u - rl -0.02 


SFG 


AGN 


12 


189.15318 


62.19889 


0.56 


11.11 


hdfn-080 


1.28c-15 


2.20e-15 


42.18 


42.42 


n 4s+ - 04 
-0.48_ 003 


-0 20+ 004 
u ' zu -0.04 


SFG 


AGN 


13 


189.27452 


62.19824 


0.90 


11.75 


hdfn-413 


<4.08e-17 


1.18e-17 


— 


40.65 


1 


-2 62+ ' 89 
z '° -0.83 


SFG 


Amb 


14 


189.45862 


62.19823 


— 


— 


— 


<4.75e-17 


<2.89e-16 


— 


— 


— 




— 


— 


15 


189.25661 


62.19619 


1.76 


12.73 


hdfn-471 


3.22c-17 


<1.09e-16 


41.79 


— 


-1 





— 


SFG 
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Table 2. - continued 



ID 


RA 


Dec. 


z 


Lm, 


HDFN ID 


/sx 


/hx 


^SX 


^HX 


HR 


/hx//r 


Spectral 


X-ray 




deg 


deg 




l°g(L Q ) 




erg/s/cm 2 


erg/s/cm 2 


log (erg/s) 


log (erg/s) 






Class 


Class 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


46 


189.03671 


62.19546 


1.34 


12.43 


— 


3.86c-17 


<3.22e-16 


41.58 


— 


-1 


— 


— 


SFG 


47 


189.05208 


62.19458 


0.28 


10.79 


— 


<5.80c-16 


<3.53c-15 


— 


— 


— 


— 


— 


— 


48 


189.22239 


62.19434 


1.27 


12.30 


hdfn-056 


9.08c-17 
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